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An iterative, CFD-based approach for aeroelastic computations in the frequency 
domain is presented. The method relies on a linearized formulation of the aeroelas¬ 
tic problem and a fixed-point iteration approach and enables the computation of the 
eigenproperties of each of the wet aeroelastic eigenmodes. Numerical experiments on 
the aeroelastic analysis and design optimization of two wing configurations illustrate 
the capability of the method for the fast and accurate aeroelastic analysis of aircraft 
configurations and its advantage over classical time-domain approaches. 


I. Introduction 

Computational aeroelasticity [T] has traditionally relied on linear methods such as the doublet 
lattice method [2] and the piston theory [3]. While these two approaches enable accurate aeroelastic 
predictions in the subsonic and supersonic regimes, respectively, they fail to produce an accurate 
description of flutter in the transonic domain. Approaches based on computational fluid dynamics 
(CFD) and Arbitrary-Lagrangian-Eulerian (ALE) formulations, have subsequently been able to 
produce very accurate flutter predictions in the transonic domain |3]. 

When used in the time-domain, CFD-based aeroelastic predictions typically rely on a post 
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processing step [5| that extracts the aeroelastic characteristics of the wet structure, such as damping 
ratio and frequencies. This post processing step can however be inaccurate and miss aeroelastic 


characteristics associated with certain modes altogether. Furthermore, depending on the initial 
conditions of the time simulation, certain modes may not be excited. 

In the present work, an alternative approach based on linearized CFD and defined in the fre¬ 
quency domain is developed. It shares some similarities with the fixed-point iteration approach 
used in the doublet lattice method that can extract flutter frequencies. However, because it relies 
on CFD, it is accurate in the transonic regime and can also extract damping ratios for non-flutter 
points. As such, it provides the capability of leading to accurate aeroelastic predictions at any 
operating point. In turn, it will be demonstrated in this paper that such a procedure can be easily 
integrated in a design optimization procedure under aeroelastic constraints. 

This paper is organized as follows. Section |TT] presents the framework of CFD-based linearized 
aeroelasticity, focusing on the frequency domain formulation. The proposed approach for the fast 
computation of the system aeroelastic characteristics is developed in Section [Till Application to the 


aeroelastic analysis and design of two wing systems is presented in Section IV Finally, conclusions 
are given in Section [V| 


II. CFD-Based Linearized Aeroelasticity in the Frequency Domain 
A. Three-field Formulation 

The following three-field ALE formulation is considered in this work to describe a fluid/structure 
interaction problem: 

• 7^(W) = 0 

' -div(cr(e(Z^)))-S =0 (1) 

-div(A =0. 

The first equation in 0 is the ALE form of the Navier-Stokes equations, the second one, the 
elastodynamics equation for the structural subsystem and the third one provides the dynamics of 
the fluid grid, represented as a pseudo-structure. All tilde notations refer to that pseudo-structure. 
In all equations, t denotes time, x{t) denotes the vector of time-dependent fluid grid position, 
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is the determinant of the Jacobian of x with respect to the reference configuration and V is the 
gradient with respect to x. 

yy denotes the vector of conservative fluid variables and U the structural displacements. T and 
TZ are the vectors of convective and diffusive fluxes, respectively, p denotes material density, e the 
strain tensor and cr the stress tensor, div is the divergence operator. S is a vector of volume forces 
acting on the structure. Finally E denotes a tensor of elasticities. 

Appropriate Dirichlet and Neumann boundary conditions as well as interface conditions apply 
to the set of coupled partial differential equations. The interface conditions between the fluid and 
structural subsystems enforce the compatibility before the two velocity fields at the interface as well 
as the equilibrium of tractions on the wet surface. Finally continuity equations are enforced at the 
wet interface between the structural and fluid mesh subsystems. The author is referred to [5] for 
more details on these interface conditions. 

B. Semi-discretization 

The fluid subsystem is discretized in space by the finite volume method and the structural and 
fluid mesh subsystems discretized by the finite element method. The fluid is assumed to be inviscid 
in the remainder of this paper. The semi-discretization, after inclusion of the appropriate boundary 

conditions, leads to the following set of coupled ordinary differential equations 

✓ 

(A(x)w)-1-(/)(w,x,x) =0 

< Mu + f“(u,u) =r’^‘(u,w) (2) 

Kx = KcU. 

X is the vector of fluid mesh nodes positions and A is a diagonal matrix containing the corresponding 
cell volumes. A dot denotes a time derivative, w is the vector of discretized conservative fluid 
variables and u the vector of structural displacements, (p denotes the numerical flux function. M is 
the mass matrix resulting from the finite element discretization and is the corresponding vector 
of internal forces, is the vector of external forces applied to the structure. Finally, the fluid 
mesh motion is modeled as being piece-wise static. K is the corresponding fictitious stiffness matrix 
and Kc is a transfer matrix representing the effect of the structural motion u at the wet interface. 
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C. Linearization 


The set of nonlinear ODEs is then linearized around an equilibrium point (wq, Xq, Uq, Uq) of the 
system. The resulting linearized system can then be written as a set of coupled linear ODEs 


AqW + Hqw + (Eo + Co)x + Gqx = 0 

Mil + Dqu + Kqu = Pqw 


(3) 


Kx = KcU. 

where the fluid mesh variable x has been eliminated through x = Su = and the linearized 

operators are 


dcj) 

Ao = A(xo), Ho = ^(wo,xo,xo), Eq = 


dA 


Cn = 


dcf) . . 

^(wo,xo,Xo) 


S, Gn = 


dcf> 

9x 


(wo,Xo,xo) 


(4) 


^£*int ^£*ext ^£ext 

Do =-^(uo,uo), Ko = -^(uo,uo)-^(wo,uo), Pq = -^^(wq,uq). 


The variables (w,u,u) denote now perturbations around the reference state (wo,Uo,Uo). For sim¬ 
plicity, the subscripts o are dropped in the remainder of this work. 


D. Structural Modal Reduction 

The dimensionality of the structural subsystem is subsequently reduced by modal truncation. 
For that purpose, a set of m eigenpairs {uJi,Zi}^^ satisfying the following generalized eigenvalue 
problem is first computed: 

Kz = w^Mz. (5) 

The eigenpairs are gathered in a diagonal matrix = diag(a;i, • • • , ujm) of eigenvalues and a modal 

basis matrix Z = [zi, • • • , z^]. 

The structural displacements and velocities are then approximated as a linear combination of 
the eigenmodes as 

u{t) « Zu„(t), v(t) « Zv,„(t). (6) 

The structural equations are then reduced by Galerkin projection as 

Imiim T H (7) 
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where fm (t) is the reduced vector of linearized aerodynamic forces defined as 


f™(t) = Z^f(t) = Z^Pw(t) = P„w(t). 


E. Frequency Domain Formulation 

The coupled equations resulting from the CFD-based ALE formulation can be written as 

liim + = PmW 

(9) 

Aw “t" Hw (E T “i- GrU^ = 0 

The aerodynamic and aeroelastic linear operators H, E, C, G and P are computed for specific 
structural configurations and free-stream boundary conditions such as free-stream Mach number 
Moo, angle of attack a as well as flight altitude h. Let s € C and assume that 


r{t) = we^*, u„(t) = 


Note that s can have a complex part. Then their derivatives are 


w(t) = swe^‘, u„(t) = su„e^‘, 


and the fluid equation becomes 


(sA + H) w + [s(E + C) + G] = 0. 


Hence 


w = - (sA + H)”^ [s(E + C) + G] u„ 


Plugging this expression for w into the structural equation leads to 


Iu„ + = -P (sA + H)"^ [s(E + C) + G] u. 


Defining the matrix of generalized aerodynamic forces A(s) = P (sA + H) ^ [s(E + C) + G], this 
leads to 


liim + + A(s)) = 0. 


Note that this equation is also valid for values of s with non-zero real part. 
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III. Aeroelastic Eigen Analysis 


A. Nonlinear Aeroelastic Eigenproblem 
Consider now a motion 


U„(t) = U/iC® 


where s = uj{i + 7 ) = iw + 017 . Then, iim = and (151 becomes 




(16) 


(17) 


This equation will now be transformed into state-space form. For that purpose, consider the 
velocity Vm{t) = Then Vh = suh- The operator .A(s) can also be decomposed into its real 

and imaginary parts as 

^(s) =^^(s)-f a^(s). ( 18 ) 

Let s be also decomposed into its real and imaginary parts s = -I- is^ (note that = W 7 and 

= w). Then 

^(s) =^^(s) + *^^(s) 




= A^{s) - 


+ s 

S-* 

Alf 


■ IS 


(19) 




= ^"(s) - 


-I- s- 


Equation 0 can then be written as 


^2„_ I I aR(\ I ('®) 


sl(su;i) + -I- A^{s) - s 


which, combined with Vh = su/j leads to the system 


sivf, 


V/i 


rj2+^«(s)- l^^^(s) 


U/1 H- r-^{sUh) = 0 


Uh + = 0 , 


= SUh 


( 20 ) 


( 21 ) 


that can be written in matrix form as 


0 

I 




U/i 

'Vh 


U/t 

V/i 


( 22 ) 


This leads to the following non-linear eigenvalue problem 


A(s) 




Uh 

Vh 


(23) 
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where A(s) is the real-valued matrix 


(s) = 






A^is) 


(24) 


Note that there are as many solutions to this eigenvalue problem as there are generalized degrees of 
freedom in the structural system. The eigenvalues come in pairs Sj = ± isj, j = 1, ■ ■ ■ ,m. 

There are two main approaches that can be followed to solve the non-linear eigenvalue prob¬ 


lem 

The first class of approaches is based on a fixed point iteration approach. In the context of the p- 
k method for solving the flutter problem, that approach was first introduced in [2] . At each iteration, 
a linear eigenvalue system is solved and the algorithm proceeds until convergence of the iterates to 
an eigenvalue. The advantage of that approach is that it relies only on a simple linear eigenvalue 
solver and is therefore straightforward to implement. A drawback of this approach is the potential 
slow convergence requiring many iterations 0. m presents an analysis of the convergence proposes 
an extension based on Newton’s method to address the slow convergence issue. That extension 
however requires the computational of the sensitivities of the aeroelastic operators with respect to 
the reduced frequency. 

The second class of approaches is based on continuation to solve the non-linear eigenvalue prob¬ 
lem. Such an approach was proposed in |8] in the context of the p-k method. The main advantage 
of this method is its robustness that guarantees a computation of all the eigenvalues. However, this 
approach is expensive as it requires computing flutter predictions for an entire trajectory from high 
altitude to the desired altitude. Further, computing the flutter path in a robust fashion may re¬ 
quire small increments that increase the computational cost. Finally, sensitivities of the aeroelastic 
operator are required for a continuation approach. 

In the remainder of this paper, the fixed point iteration approach is chosen for its simplicity of 
implementation. It is described in Section [ilIB[ 
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B. Fixed-point Iteration Algorithm 


The algorithm proceeds, for each of the m structural eigenvalues, by iteratively computing one 


solution of the eigenvalue problem (23). For that purpose, an initial gness is provided by the dry 
eigenvalues and associated eigenvectors. The operator A is then evalnated at the dry mode and 
all of its eigenvalues computed. The eigenvalue are sorted according to their imaginary parts and 
the eigenvalue associated with the index of interest retained. The operator A is then evalnated 
for that eigenvalue. The procedure is repeated until convergence. A pseudo-code is provided in 
Algorithm If convergence is not achieved after a maximnm number of allowed iterations, the 
continuation procedure mentioned above is triggered. 

The main difference with the p-k method developed in [2] is the fact that the real part of the 
eigenvalues is retained during the iterations, leading to a possible convergence of the procedure to 
an aeroelastic eigenvalne with non-zero real part. This means that aeroelastic eigenvalue can be 
computed for non flutter configurations as well. This is due to the fact that the linearized CFD- 
based framework developed in Section is equally valid at flutter and non-flntter conflgnrations, 
unlike the p-k method that is based on the donblet lattice method, valid only for harmonic motions. 


Once the eigenvalues {sj = nonlinear eigenproblem (23) are computed, the 


associated aeroelastic frequencies and damping ratios can be readily computed as 

^ T *1 . 1 

Vj = - ; „ , J = 1, - • • .W. 


O’ 




(25) 


C. Complexity analysis 

The main cost associated with the proposed procednre originates from the computation of 
the linear operator A(s) for values of s arising in the procedure. Computing A(s) reqnires the 
computation of the matrix of generalized forces A(s) = P (sA + H)~^ [s(E -f C) -|- G]. In practice, 
computing the term (sA -|- H)~^ [s(E -t- C) -I- G] is the most expensive step. It relies on the solution 
of m large-scale, sparse systems of equations (sA -t- H) = [s(E -t- C) -f G] e^, i = 1, • • • , m. These 
solutions are typically computed by a Krylov-subspace iterative technique. In the present work, the 
pre-conditioned Generalized Minimal Residual technique (GMRES) [HI is used. 

Assuming, for simplicity, that the computation of each of the m aeroelastic eigenvalues requires 






Algorithm 1 Iterative computation of the eigenvalues of the nonlinear eigenproblem (231 


Require: Dry structural frequencies tolerance e 

Ensure: Eigenvalues {sj = sf + 

1: for j = 1, - ■ ■ ,m do 
2: Let Sj = iojj and A: = 0 

3: while fc = 0 or sj — > e do 

4: A: = fc + 1 

5: Compute all the eigenvalues {A; = Xf ± iXj}JLi of A(s*“^) and order them by frequency such that 


A{ < ■ • ■ < A^ 


5; Let — A^ iXj 

6: if A: > Ajmax then 

7: Call the continuation procedure 

8: end if 

9: end while 

10: end for 


Niter iterations, Niterm? such solutions of linear systems are required. 


D. Computation of the Flutter Speed Index 

The Flutter Speed Index (FSI) is a non-dimensional quantity that characterizes the flutter 
properties of an aeroelastic system. It is computed for a fixed free-stream Mach number Mao and 
density pao- The free-stream pressure Pao is increased until the onset of flutter. This is equivalent, 
for a perfect gas, to increasing the free-stream velocity Vao as 


14o = 


7^00 


Ma 


(26) 


with 7 = 1.4. The FSI is then defined for the critical value of flutter velocity 1/^'' as 


FSI = 


yci- 
* oo 


(27) 


bsUJar/Ji 

where bs is a reference length such as the semi-chord of the wing at the root, Wq, is the frequency 
of the first torsional dry mode of the structural system. The parameter p, is defined as 


ms 



(28) 
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where ms is the structural mass of the system and V is the volume of a conical frustum with lower 
and upper based diameters being the stream-wise chords at the root and the tip of the wing and 
with height equal to the semi-span of the wing of the system. 

Flutter can be in practice detected from the outputs of Algorithm by the presence of an 
eigenvalue Sj of the operator A with positive real part. In practice, Algorithmis run for increasing 
values of the free-stream pressure until flutter is detected. 

IV. Applications 

A. Flutter Analysis of the Agard Wing 445.6 

1. Agard Wing 445-6 model 

The weakened Agard 445.6 model 3 research wing model is described in m- The dimensions of 
the wing are reported in Figure A finite element model (FEM) consisting of 800 shell elements, 
corresponding to 2646 degrees of freedom (dofs) is created and the eigenfrequencies associated with 
its first m = 4 structural modes compared in Table to the ones measured experimentally in |10| . 
Good agreements are observed between the two sets of eigenfrequencies. 

A finite volume mesh with more than 22,000 nodes is then created, leading to more than 110, 000 
dofs in the CFD model. 


Mode Frequency (Hz) Frequency (Hz) Relative 

number (experimental model in [10] ') (present FEM model) discrepancy 


1 

9.60 

9.84 

2.5% 

2 

38.17 

39.54 

3.6% 

3 

48.35 

50.50 

4.4% 

4 

91.54 

96.95 

5.9% 


Table 1 AGARD Wing 445.6 retained modes for the structure 


2. Flutter Analysis 

A flutter analysis is first carried on for three different Mach numbers M^o G {0.8,0.95,1.1} and 
three altitudes h G {0,10000, 20000} ft. For every combination of flight conditions, the proposed 
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Fig. 1 Geometry of the Agard wing 


iterative approach is applied to predict the aeroelastic eigen-characteristics of the system. Hence, 
for each case, four sets of real and imaginary parts for the eigenvalues are computed, as well as the 


corresponding damping ratio, as defined in Eq. (25l. The densities and pressures associated with 


each altitude are given in TableThe results are graphically depicted in Figure]^ a) for aeroelastic 
predictions at sea level, in Figure [^b) for h = 10,000 ft and in Figure [|;c) for h = 20,000 ft. 
In addition, the results are compared to predictions in the time domain for which damping and 
frequencies characteristics are extracted by fitting a canonical response to the transient lift. This 
simple approach only allows the extraction of a single eigenvalue, as opposed to the approach 
proposed in this paper that allows the extraction of all the eigenvalues. The comparisons with 
the proposed iterative approach show that the time domain are able to accurately predict the 
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characteristics of one of the aeroelastic eigenmodes, but that eigenmode is not always associated 
with the smallest damping ratio. As a result, the time domain approach ofter overestimates the 
minimal damping ratio of the system, as observed in Figurej^where the damping ratios computed by 
each approach are reported. Furthermore, the time-domain mispredicts flutter for M^o = 0.95 and 
h = 10k ft. On the other hand, the proposed iterative approach returns all the eigenvalues and hence 
the one associated with smallest damping. The aeroelastic analysis for Moo = 0.95 and h = 20, 000 
ft takes 4.25 minutes on 8 processors for a tolerance e = 10“®. Computing the steady-state takes 
1.1 minute and the iterative approach takes 3.15 minutes. 


Altitude (ft) Pressure (Ibf.in^) 

Density (Ibf.s^.in ‘^) 

0 (sea level) 

14.7 

1.1463 X 10"'^ 

10,000 

10.11 

0.8466 X 10"'^ 

20, 000 

6.76 

0.6112 X 10"'^ 


Table 2 Free-stream pressures and densities for the altitudes considered 

Next, FSIs are computed for a range of Mach numbers varying from Moo = -499 to Moo = 1.141 
using the proposed approach. For that purpose, the free-stream density is first fixed to its sea level 
value. The results are reported in Figure One can observe a characteristic flutter dip around 
Moo = 0.96 and a range of flutter conditions for Moo G [0.87,1.05]. 

B. Analysis and Design Optimization Based on the ARW-2 Wing 

1. Aeroelastic Research Wing ARW-2 Model 

The second model is a more realistic wing. It is the Aerodynamic Research Wing Number 2 
(ARW-2) which was designed for flutter suppression and gust load alleviation [TTl IT^. The wing 
is based on a supercritical airfoil, has an aspect ratio of 10.3 and a leading-edge sweepback angle 
of 28.8°. The geometry of the wing is depicted in Figure A FE model of the wing consists of 
1765 shell, beam and truss elements, resulting in 2736 degrees of freedom. The structural FE model 
is depicted in Figure m = 8 structural eigenmodes are then selected to model the structural 
subsystem. The corresponding eigenfrequencies are reported in Table 
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Fig. 2 Comparison of the aeroelastic eigenvalues predicted by the frequency and time domain 
methods for the Agard wing 


The fluid computational domain consists of an unstructured grid with 47, 847 vertices corre¬ 
sponding to more than 260,000 dofs. Part of the fluid mesh is shown in Figure 


2. Flutter Analysis 

A flutter analysis is carried on for the three Mach numbers Moo G {0.8,0.95,1.1} and altitudes 
h S {0,10000, 20000} ft. For each combination of flight conditions, the proposed approach is applied 
to predict the aeroelastic eigen-characteristics of the system. 

Hence, for each case, eight sets of real and imaginary parts for the eigenvalues are computed. 
The eigenvalues are reported in Figure [^a) for aeroelastic predictions at sea level, in Figure |^b) for 
h = 10,000 ft and in Figure [^c) for h = 20,000 ft, together with predictions in the time domain. 
One can observe that the time domain approach is again associated with large approximation errors 
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Fig. 3 Comparison of the smallest aeroelastic damping ratios predicted by the frequency and 
time domain methods for the Agard wing 


of the aeroelastic eigenvalues, especially at M^o = 0.95 and Moo = l-l- In turn, the minimal 
damping ratio is often mispredicted by this time-domain approach,as observed in Figure The 
aeroelastic analysis for Moo = 0.95 at sea level takes approximatively 6.8 minutes on 32 processors for 
a tolerance e = 10“®. Computing the steady-state takes 1.2 minute and the iterative eigenanalysis 
approach takes 5.6 minutes. 

Finally, the flutter speed indices associated with the wing at sea level altitude are computed and 
reported in Figure |10| One can observe that flutter is never reached, even at sea level, confirming 
the proper flutter-free design of the ARW-2 wing. 
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Fig. 4 Variations with the Mach number of the Flutter Speed Index and comparison with the 
sea level conditions for the Agard wing 



Fig. 5 Geometry of the ARW-2 wing 


3. Design Optimization 

Next, the proposed aeroelastic analysis approach is included in a design optimization procedure 
in order to design an aircraft system under specific dynamic and static aeroelastic constraints. For 
that purpose, the ARW-2 wing is considered as a baseline design at M^o = 0.95 and sea level 
altitude. In this study the first m = 6 structural modes are considered.The following optimization 
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Mode number Frequency (Hz) 

1 

7.92 

2 

28.83 

3 

33.27 

4 

57.04 

5 

73.70 

6 

86.08 

7 

93.91 

8 

110.87 


Table 3 ARW2 retained modes for the structure 



Fig. 6 Finite element model of the ARW-2 wing 



(a) Surface and symmetry plane 


(b) Fluid surface mesh 


Fig. 7 ARW-2 wing fluid mesh 
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Fig. 8 Comparison of the aeroelastic eigenvalues predicted by the frequency and time domain 
methods for the ARW-2 wing 


problem is considered 


L{fi) 
max —;—- 

fiGV D(fi) 


S.t Ty(/r)<Wma: 
^(m) — ^max 


(29) 


where fi G denotes a vector of design variables and T> C is the design space. and 

Z?(/r) respectively denote the lift and drag associated with the static aeroelastic equilibrium for the 
design configuration /x. denotes the weight associated with the design /x. The vector cr(/x) 

contains the nodal values of the von Mises stress in the structure at that configuration. Finally, 
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Fig. 9 Comparison of the smallest aeroelastic damping ratios predicted by the frequency and 
time domain methods for the ARW-2 wing 

ri{iJ,) is the vector of m modal damping ratios associated with the configuration for the flight 
condition of interest. Wmax) o-max and provide respective bounds on the weight, stress and 
damping ratios. 

Tthere are = 6 design variables. The first three are associated with the shape of the wing 
and respectively control the sweep, twist and dihedral angles. The last three variables control the 
thickness of three distinct subsets of spar and rib elements in the wing. The design space V is a 
box domain defined in Table Note that all six design variables are non dimensional and scaled. 

The nonlinear program ( [2^ is solved with Matlab’s fmincon active set method with convergence 
tolerances set to 10“®. The gradients are computed analytically except for the gradient of the 
damping ratio that is computed by finite differences. This requires a total of + 1 = 7 calls to 
the iterative procedure introduced in this paper for each design considered. 
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Fig. 10 Variations with the Mach number of the Flutter Speed Index and comparison with 
the sea level conditions for the ARW-2 wing 


Design variable 

Nature 

Lower bound Upper bound 

Si 

Sweep angle 

-0.1 

0.1 

S2 

Twist angle 

-0.1 

0.1 

S3 

Dihedral angle 

-0.1 

0.1 

S4 

Element thickness 

-0.1 

0.1 

Ss 

Element thickness 

-0.1 

0.1 

Sq 

Element thickness 

-0.1 

0.1 


Table 4 Design domain for the ARW-2 wing 


The initial design is the one associated with the original ARW-2 wing. Its physical properties 
are reported in Table The bounds in (291 are set to Wmax = 400 lbs, CTmax = 2.5 x 10“^ and 
??min = 2 X 10“^. As such, the dynamic aeroelastic constraint is violated for the initial design. The 
optimized design is found after 13 major iterations and reported in Table and the two design 
configurations depicted in Figure |11[ Unlike the original design, the damping ratio constraint is 
satisfied and active for the optimized design. The lift/drag ratio has increased as well as the weight. 
The maximum von Mises stress has decreased. 

The overall design procedure takes 8 hours on 32 CPUs. About 1 hour is associated with the 
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Property Initial design Optimized design 


Si 

0 

-0.1 

S2 

0 

-0.1 

S3 

0 

0.1 

S4 

0 

0.1 

ss 

0 

0.1 

S6 

0 

0.011 

L(/x) 

1.283 

1.355 

D(p) 



W{fi) (lbs) 

349.9 

365.4 

max(cr(/x)) (psi) 

1.54 X 10'‘ 

1.28 X 10' 


1.82 X 10-^ 

2.00 X 10 


Table 5 Properties of the initial and optimized designs 


computation of the multiple static constraint and its derivatives and 7 hours to compute the dynamic 
constraints and their sensitivities. 


The variation of the design over the optimization iterations is reported in Figure 12 and the 
corresponding objective function in Figure |13| 


V. Conclusions 

This paper presents a practical and accurate approach for determining the aeroelastic charac¬ 
teristic of an aircraft configuration. The approach relies on a linearized CFD-based formulation 
together with fixed-point iteration procedure. The application of the proposed approach to the 
aeroelastic analysis and design optimization of two realistic wing configurations illustrates its capa¬ 
bility to enables fast and accurate aeroelastic predictions in the subsonic, transonic and supersonic 
domains. 
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Fig. 12 Variations with the optimization iteration of the design variables 


pp. 856-869. 

[10] Yates, E. C., “AGARD Standard Aeroelastic Configurations For Dynamic Response - 1 - Wing 445.6,” 
NASA, 1987. 

[11] Seidel, D., Sandford, M., and Eckstrom, C., “Measured unsteady transonic aerodynamic characteristics 
of an elastic supercritical wing,” Journal of aircraft, Vol. 24, No. 4, April 1987, pp. 225-230. 

[12] Sandford, M., Seidel, D., and Eckstrom, C., “Geometrical and structural properties of an aeroelastic 
research wing (ARW-2),” NASA Technical Memorandum, May 2003, pp. 1-39. 


22 







